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A SYSTEM FOR COMPACTION AND RECONSTRUCTION OF WAVELET DATA 

HELD OF THE INVENTION 

TTie piesent invention lelaies lo data processing systenos fenerally and to 
systems utilizing wavelet transfoims in particular. 

BACKGROUND OF THE INVENTION 

Data processing typically involves tiansfonning the data &om one space to 
another. The best known transfonnaticai is the Fourier transfonnati(»i which transforms a 
signal from a time space to a frequency space, indicating thereby the frequency components 
of the input signal 

Recently, wavelet transforms have been developed where each wavelet 
transform has a dilferent curving "wavelet" shape and the result of the transformation 
indicates to what extent the input signal, a digital signal, is formed of a combination of the 
different wavelets. The article "Wavelets and Signal Processing", by Olivier Rioul and Martin 
Veiterli, published in IEEE SP' Magazine. October 1991, pp. 14 - 38, describes prior art 
wavelet transform operations and is incorporated herdn by reference. 

Typically; wavelet transforms are implemented in a multi-resohition manner. 
The prior art mulli-resolution method and exemplary output, for a one-dimensional signal, are 
illustrated in Figs. 1 A and 2, respectively, to which reference is now made. 

Shown in fig. lA are three "wavelet transform units" 10a - lOc. Each wavelet 
uansform unit lOa - 10c is formed of its respective wavelet low pass filler (WLPF) 12a - 12c, 
wavelet high pass filter (WHPI') 14a - J4c and samplers 16a - 16c. The wavelet low and 
high pass filters 12a - 12c and 14a - 14c break the signal into its high and low wavelet 
frequency components and the samplers 16a - 16c sample the output of each filtei, thereby 
changing the resolution of the signal For one-dimensional input signals, the sampleis 16a - 
16c typically take every second datapoint, thereby reducing the number of points by two with 
each level Other sampling rates can also be implemented. 

Each transform unit lOa - lOc operates on a different resolution signal. 
Therefore, the elements of each transform unit are labeled to indicate the resolution level or 
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iteiation. 1 - 3. at which they operate, wherein the highest level has the lowest resolution 
(least number of points). The first transform unit 10a operates on the input signal and the 
other transform units 10b - 10c operate on the low frequency output signals of the transform 
nnil 10a - 10b. respectively, of the level previous to it The signals of Kg. 1 are labeled to 
5 indicate the level and whether they are the low (L) or high (H) frequency output of the level 
The output of a wavelet transform operation comprises the following signals, 
the high frequency output signals HI - H3 and the low frequency output signal L3 of the last 
level The low frequency output signal of the highest level, in this example signal L3, is 

herein called the 'core" signal. 
,0 An exemplary input signal, labeled 18, and its conesponding three high 

frequency output signals are shown in Fig. 2. As can be seen, the resolution is reduced at 
each level Furthermore, the activity of the lower resolution signals is near at the inflection 

point of the input signal 18. 

If the input signal 18 is an image or some other two-dimensional signal, the 
15 wavelet low and high pass filters can either be separable or non-separable. The operations 
for separable and non-separable filters are shown in Figs. IB and IC, respectively, to which 

reference is now made. 

The wavelet transform units 10 for separable filtering include low and high 
pass wavelet filters 11a and lib (Fig. IB) and samplers 13a for rows and low and high pass 

20 wavelet filters 1 Ic and 1 Id and samplers 13b for columns. The rows of the input signal first 
pass through the low and high pass wavelet filters 11a and lib and are sampled by the 
samplers 13a, thereby reducing the number of pixels in each row by two. The columns of 
the row filtered signals then pass through the low and high pass wavelet filters 11c and lid 
and are sampled by the samplers l3b, thereby reducing the number of pixels in the columns 

25 by two. Alternatively, the columns can be operated on first, followed by the rows. 

The separable filtering method produces three high frequency signals Ha, Hb, 
He and one low frequency signal L at each level Thus, the first level produces the signals 
Hla, Hlb, Hlc and LI. Each signal has one-fourth the number of pixels as the input signal 
The low frequency signal LI is then provided as the input signal to the next wavelet 

30 transform unit 10. 

The wavelet transform units 10 for non-separable filtering operate differentiy 
for even- and odd-numbered levels, as illusuated in Rg. IC. For odd levels, the input signal 



is operated on by non-sepaiaWe low and high pass wavelet filteis 15a and 15b, lespectively, 
and quincunx odd sampleis 17a. Odd sampleis 17a perform "quincuux" subsampling which 
removes every even datapoint in odd numbered rows and every odd datapoint in even 
numbered rows, leaving a checkerboard-like output image, illustrated in part in Fig. lOR In 
5 order to mainlab the image, the datapoints are not actually removed; instead, they just arc 
not operated upon by the foDowing even level 

For even levels, the wavelet transform unit 10 operates on the checkerboard- 
like image with non-separable low and high pass wavelet filters 15c and 15d, respectively, 
and quincunx even samplers 17b. Even samplers 17b also perform quincunx sampling on the 
10 checkerboard-like image, resulting in a more noraial looldng image, shown in part in Fig. 
lOA, which is the input to the odd-level wavelet transform unit The image produced by even 
samplers 17b has one-half the number of elements as the input image to Ihe odd samplers 
17a. 

The anicle "Predictive Inlerscalc Image Coding Using Vector Quantization", 
15 by M. Antonini el al and published in Sipal Process ing V:Theories and Applications., L. 
Tones, E Masgrau and M.A. Lagunas, editors, Elsevier Science Publishers, B.V. 1990. pp. 
1091 - 1094, describes an image compression method using biorthogonal wavelet transforms 
and an interscale prediction scheme. The prediction scheme predicts the location and 
amplitude of edges in an image at one resoluUon based on their locations and amplitudes at 

20 the next highest resolution. 

The article "Image Coding using Lattice Vector Quantization of Wavelet 
Coefficients", by M. Antonini et al., ICASSP 91, pp. 2273 - 2276 describes a compression 
method for wavelet transform data. 

U.S. Patent 5,014,134 to Lawton et al. describes an image compression and 
25 decompression system which utilizes wavelet transforms. 
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SUMMARY OF THE PRESENT INVENTION 

II is an object of the present invention to provide a system and method of 
compaction and icconsiniction of wavelet tiansfonn data which preserves areas of high 

5 iBformaiion coBtcnt in an inpait signal. The input signal can be a one-dimensional signal, a 
two-dimensional image or a many dimensional signal. 

There is therefore provided, in accordance with a prefened embodiment of the 
present invention, a system for compaction and reconstruction of wavelet transform data of 
an input signal where the wavelet transform data is formed of a pluraUty of resolution levels. 

10 The system includes a compaction unit and a convener. The compaction unit is operative on 
forests of trees, wherein each Uee has branches and each branch relates one datapoint in a low 
resolution level with the entirety of its descendant datapoinis in higher resolution levels. The 
compaction unit prunes branches from the trees and trees from the forest to remove aH but 
areas of high information content in the input signal, thereby to produce pruned trees. The 

15 converter converts the pruned uees to reconstructed wavelet transform data. 

Additionally, in accordance with a prefened embodiment of the present 
invention, the compaction unit includes a) a non-null pruning unit for removing branches of 
the trees having datapoints thereon whose wavelet uansform values are below predetermined 
threshold levels, producing thereby non-null pruned uees and b) a zero-crossed pruning unit 

20 for determining locations of transients in the input signal and for creating zero-aossed pruned 
trees having therein only datapoints of the non-null pruned uees which are at or near the 

locations of uansients. 

Moreover, in accordance v^th a preferred embodiment of the present invention, 
the datapoints of the zero-crossed pruned uees form nodes and the compaction unit further 
25 includes apparatus for associating datapoints of the non-nuD pruned uees with nodes of said 
zero-crossed pruned uees, wherein each datapoint can be within a local area of many nodes 
but is associated with only one of them. 

Further, in accordance with a prefened embodiment of the present invention, 
the system includes a nullifier, operative prior to the compaction unit, for nulling those 
30 datapoinis of the wavelet uansfonn data which are below the predeiennined threshold levels. 

Stin further, in accordance with a prefened embodiment of the present 
invention, the zero-crossed pruning unit includes apparatus for detennining locations of zero- 
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crossings in the wavelet tiansfonn data. 

Finally, the method peifonns the operations of the elements of the system. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



The present invention will be understood and appreciated more fully from the 
foDowing detailed description taJcen in conjunction with the drawings in which: 
5 Fig. lA is a block diagram illustration of a prior multi-resolution wavelet 

transform operation; 

Fig. IB is a block diagram illustration of one section of a prior art wavelet 

transform operation for separable filtering; 

Fig. IC is a block diagram illustration of two sections of a prior art wavelet 

10 transform operation for non-separable fUtering; 

Fig. 2 is a graphical illustration of an input signal and versions of the signal 
after wavelet uansformation having high frequency information therein; 

Fig. 3 is a block diagram illustration of a wavelet transform compaction and 
reconstruction system, constructed and operaUve in accordance with a preferred embodiment 

15 of the present invention; 

Fig. 4 is a graphical illusuation of the parent-child relationships between the 
signals of Fig. 2, useful in understanding the operation of the present invention; 

Figs. 5A and 5B are schematic illustrations of fuU and pruned trees, useful in 

UDdersianding the operation of the system of Fig. 3; 
20 Fig. 6 is a block diagram illustiaiion of a method of nullifying some wavelet 

transform data, useful in the system of Hg. 3; 

Fig. 7 is a flow chart illustration of the operation of a compaction unit forming 

part of the system of Fig. 3; 

Fig. 8 is a flow chart illustration of a method of producing a zero-crossed tree, 

25 forming part of the operations of the compaction unit of Fig. 7; 

Fig. 9 is a flow chart illustration of a method of computing zero-crossing 
locations forming part of the method of Fig. 8; 

Fig. lOA is a schematic illustration of a local area within a two-dimensional 

image; 

30 Fig. lOB is a schematic illustration of a local area for an odd level wherein 

quincunx subsampling was performed; and 

Fig. 1 1 is a flow chart illustration of a method of producing a fully pruned tree 



foiming part of the method of Hg. 7. 



DETAILED DESCRffTlON OF A PREFERRED EMBODIMENT 

The present inveniion will be described hercinbelow in the context of one and 
5 iwo-dimensional input signals, it being understood that input signals of other orders of 
magnitude aie also included within the scope of the present invention. 

Reference is now made to Fig. 3 which illustrates a wavelet uansform 
compaction and reconstruction system, constructed and opcraUve in accordance with a 
prefened embodiment of the present invention. The system of the present invention typically 
10 comprises a compaction system 20 for compacting wavelet transform signals of a digital input 
signal and a reconstruction system 22 for reconstructing the input signal from the compacted 

wavelet transform signals. 

The compaciion system 20 typically comprises a wavelet transformer 24, a 
nulEfier 26 for zeroing out values of the wavelet transform data W which are below a 
15 predetermined human-defined threshold and a compaction unit 28 which compacts the output 

W*' of the nullifier 26. 

The vwivelet transformer 24 comprises a predetermined plurality of wavelet 

transform units 10 combined as described in the Background of the Present Invention. For 

most of the examples herein, wavelet transformer 24 comprises three wavelet uansform units 
20 10, however it will be understood that other numbers of miits 10, typically in the range of 

three to eight, are suitable. 

The wavelet transform units 10 can be separable or non-separable. It is noted 

that the type of wavelet transform utilized (either separable or non-separable) affects the 

percentage of compaction to be achieved, as will be described hereinbelow. 
25 As is known in the art, the wavelet transformer 24 produces two types of 

signals, the core signal and the detaDed signal. The core sipal is the low frequency signal, 

'labeled L3 in Fig. 1, which has a low resolution. The detaflcd signal comprises the high 

frequency signals HI - H3 each having a different resolution. 

If separable filtering was utilized on a two-dimensional input image, each high 
30 frequency signal HI - H3 is formed of three signals. For the remainder of this discussion, 

the term HI will be utilized and it will refer to the output of both separable and non-separable 

filtering unless otherwise specified. 
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•n,e compaction miil 28 l^picaDy deiennine the aicas of "Tugh infonnation 
content- oi of detail in the nnllified wavelet txansfonn data W«. H.e«e areas can be defined 
as those areas which have transients in them and which exist in most or aU of the high 
frequency signals of the nullified wavelet data W". 

b accordance v^th the present invention and as illustrated in Figs. 4. 5A and 
5B to which reference is now briefly made, the compaction unit 28 organizes the detailed 
signal as a "forest" of "trees" com.ecting the datapoints of the higher resolution signals (HI 
aBd H2) with their corresponding datapoints in the lowest resolution signal (H3). TTie 
datapoints of the lowest resolution level form the "roots" of the trees and therefore, the lowest 
resolution signal H3 is herein denoted the "root signal". There is no more than one tree for 

each datapoint in the root signal 

The compaction unit 28 then "prunes" the uees to remove datapoints having 
UtUe information content The compaction unit 28 can also remove trees from the forest if 

there is no information in a tree. 

Figs. 5A and 5B illustrate two full Uees 30 and two pruned trees 32 for two 
datapoints (0 and (i+ 1) of the root signal of a one-dimensional signal. In this example, there 
are five levels and thus, five high frequency signals. HI - H5. are illustrated. In this example, 
each wavelet tiansfonn unit 10 reduced the resolution of the signal by % Therefore, each 
daupoint at one resolution level corresponds to two datapoints in the next highest resolution 

20 level 

It is noted that each tree 30 and 32 comprises a multiplicity of branches 34 
connecting each datapoinl or node 35 of a signal with its conesponding daupoints of the 
signal having the next highest resolution. It is also noted that there are no branches 34 of the 
pruned trees 32 which are not comiectcd to the structure. Thus, analogous to pruning a 
25 physical tree, when one branch 34 is removed, aD of the branches 34 connected to it in the 
higher resolution signals must also be removed. 

Referring back to Fig. 3. the compaction unit 28 produces a forest of pruned 
trees 32. at most one tree for each datapoint in the root signal. The forest of pruned trees 32 
contains' generaDy all the essential information in the input signal, but in a more compact 
format than the output of the wavelet uansfonn. The compact forest of pruned Uees 32 can 
then be uansmitted as desired and as illustrated in Fig. 3. For the purposes of the present 
discussion, "uansmission" indicates transmission between any two elements, whether through 
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commimication devices or by storing on a storage device. 

He leconstrucUon system 22 of the piesent invention typically comprises a 
convener 40 for converting the pruned forest hack to a wavelet transform representation 
and an inverse wavelet transformer 42 for performing the inverse wavelet transform operation, 
thereby to reconstruct the input signal from the wavelet transform representation W. 
Because the compaction system 20 maintains generally the essential details of the input signal, 
the reconstructed signal suongly resembles the input signal. 

It will be appreciated that the system of the piesent invention can be utilized 
iB any situation where the input signal must be reduced in size while the information content 
v^thin it be maintained. For example, if the input signal is a high quality digital image, 
comprised of more than one-quarter million picture elements or pixels, it can be compacted 
with the present system for storage and reconstructed for operation. Alternatively, it can be 
transmitted to a second location after compaction and it can be reconstructed at the second 
location. 

The compacted form of the wavelet transform lends itself easfly to further 
compression, as foDows: The core signal may be further compressed by known Discrete 
Cosine Transform (DCT) methods; the ftiD-pruned forest topology may be encoded by 
standard tree encoding methods; and the datapoints of each node of the fully-pruned forest 
may be considered as a vector and thus, vector quantization may be applied to compress 
them. 

Reference is now made to Fig. 6 which illustrates the operations of nullifxer 
26. The equations performed by the nuUifier 26, for a two dimensional input signal (e.g. an 
image) are listed hereinbelow in the section labeled "NuUifier Equations". The index « 
indicates the resolution level and (ij) are spatial indices. The nullifier 26 determines a 
threshold level for nullilybg. based on the type of input signal and on its particular features. 
Any values of the wavelet data less than the threshold value are nullified (e.g. set to 0). 

Because the wavelet transform data W is mulli-resolutional and the values in 
a signal at one resolution level are not equivalent to those in a signal at another resolution 
leve^ due to human response mechanisms to the specific type of input signal (audio, video, 
etc), the data W, in step 52, has to be normalized according to human responses. 

The human response function R(f depends on the resolution level { and also 
on the position of each datapoint (ij) within the image. For the type of wavelet transforms 
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described beieinabove, an exemplary human response function R(f,iJ) i 
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R(«.iJ) - m - 0.95*2.2-« 

Prior to noroializaiion, in step 50. the absolute values of the data arc 
determined, producing signal W,. The output of step 52 is a signal W,. 

After normalization, the data is quantized (step 54). vvith a quantization level 
Q. to reduce the number of possible values, and then clipped (step 56) to maintain the number 
of possible values to a finite and reasonable amount, such as 1024. 

It wiU be appreciated that steps 50 - 56 can be performed at once, if desired. 
A histogram of the resultant normalized, quantized data from all of the high 
frequency signals, is produced in step 58. b step 60 an overall threshold level T. is 
determined, based on the desired percentage, such as 80%. of the histogram to be nullified. 
The overaU threshold level is the value of the histogram for which the desired percentage of 
15 the data has values less than or equal to it 

The overaU threshold level cannot be utilized for the threshold operation (step 
64) since its value is a normalized one. mrefore. in step 62. the threshold is -denormalized" 
(i.e. the inverse of the normalization operation) to provide threshold values T(«,iJ) for each 

resolution level t and across the image. 

In step 64, the threshold operation is performed wherein each datapoint in each 
high frequency signal HI - H3 of the wavelet transform data W is compared to its 
conesponding threshold and those datapoints uWch are less in absolute value than the 
corresponding threshold are nullified. 

Reference is now made to Hg. 7 wbch illustiates the operations of compaction 
25 unit 28. Compaction unit 28 performs two types of pruning, a first pnming (step 70) to 
remove those branches which have nuU date on them, and a second pruning (step 72) to 
remove those branches which do not have transient date. Since transients are typicaDy more 
than one datepoint wide, compaction unit 28 then enlarges (step 74) the nodes of the resultant 
pruned tree to include at least some of the non-null datapoints near the datepoint of the node. 
30 The pruning steps 70 - 74 can be performed by creating lists of the nodes 

which belong in each pruned tree. In accordance with an altemaave embodiment of the 
present invenUon, each datepoint can have many attributes and pruning is performed by 

11 
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updating the attributes of each node which belongs in the pruned tiee. The steps 70 - 74 will 
be described as though new pruned uees aie produced, it being understood that either method 
can be performed. 

Prior to discussing the specific details of each operation of compaction unit 28, 
5 ^ define a parent node and its child nodes of the level below it For a one-dimensional 
signal, a parent node W»(Ci) has the following children nodes: ^(t-l,2i^l), V^(t-\,2i). 

For a two-dimensional signal upon which separable filtering was performed, 
a parent node W°(J.d.ij) has the foDowing child nodes: W«(M.d.2i-^1.2j+l), ^Hi- 
l.d,2i-H,2j). W<'(M,d.2i.2j+l) and V^{t-WU2!). where d indicates the type of filtering 
10 plr'form^. ' For example, d - 1 indicates that the lo^ were low pass filtered and the 
columns were high pass filtered, d - 2 indicates that both rows and columns were high pass 
mtered and d = 3 indicates that the rows were high pass filtered and the columns were low 
pass Altered. 

For a two-dimensional signal upon which quincunx subsampling was 
15 performed, a parent node W«(f.iJ) of an even-numbered level has the following child nodes 
of the odd-numbered level below it: WO(M.2i.2j). W<'(M.2i-H,2j+l). A parent node W«(tij) 
of an odd level has the foUowing child nodes of the even-numbered level below it: W(M,i j) 
and W°(M,i+l j) for j even and W«(«-l.ij) and W«(M,i-lj) for j odd. It is noted that in the 
odd levels, if i+j is odd, the datapoint docs not exist 
20 Similar relationships can be built for input signals of higher dimensions. 

Ib step 70, the compaction unit 28 reviews the entirety of nodes in the forest 
of trees. For each node which is null and whose descendants do not have significant 
information in them, where "significance" is defined beieinbelow in equation (2). the node 
and its descendants are removed from the uee to which it belongs. TypicaDy, the compaction 
25 unit 28 begins the review with the root signal H3. 

The significance S of the data in the descendants of the node is determined as 

follows: 
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S - 23 W«(*)R(*) 

where (*) indicates indices and denotes a descendant of the node 35. W«(*) is the value of 
the descendant, R(*) is the value of the noiroalization function for the descendant, and the 
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sum 



is over all the descendants of the node. S has to be above a piedeiennined threshold 
value in order for the node to be kept within the pruned tree stnictiue. 

For those nodes 35 which have nuD values, the compaction unit 28 determines 
whether or not the branches 34 attached to the node have significant data in them. If so. it 
saves the node as part of the pruned tree structure. 

The compaction unit 28 repeats the above determination for each of the 
datapoints in the root signal and for each of the high frequency signals. The result is a forest 

of "non-null* pruned trees. 

In step 72. the compaction unit 28 reviews the non-null pruned tree and the 
wavelet transformation W to determine the areas conesponding to "uansicnts". A transient 
is herein defined as the location where the wavelet transformation value changes from positWc 
to negaUve or vise-versa; in other words, a transient is the location of a -zero-aossing*. By 
choosing the appropriate wavelet bases for the wavelet uansformation operation (performed 
by wavelet transformer 24 (Fig. 3)), the 2cro crossings can be forced to be at locations of 
maximum transitions in the input signal 

Selection of appropriate wavelet bases is discussed in the foDowing article 
which is incorporated herein by reference: Stephane Mallet, "Zero-Crossings of a Wavelet 
Transfonn", TFFF Transactipt»<: nt. Tuformation Theory Vol. 37, No. 4, July 1991, pp. 1020 - 
1033. 

The operations performed are provided in the flow charts of Figs. 8 and 9, to 
which refeience is now made. Fig. 8 illusirates the operations performed to produce the zero- 
crossed pruned trees and Fig. 9 illustrates the operations performed to determine where the 

zero-crossings occur. 

In step 80 of Fig. 8, which is deiaUed in Fig. 9, the compaction unit 28 
determines which nodes 35 of the non.-null pruned forest are part of transients. The locadon 
of the transients are determined through an iterative process operating on the non-nuUified 

wavelet uansform W. 

The iteration process shown in Fig. 9 operates on a "local area" or 
neighborhood in W of each node 35 in the non-null pruned tree. Local areas for a two- 
dimensional input signal are shovm in Figs. lOA and lOB to which reference is now briefly 
made. For signals upon which two-dimensional separable sampling was performed, the 
neighborhood is typically an NxN square of data surrounding the cunent node 35, where N 
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is typicaDy 3 - 7. This is illusirated in Fig. lOA. If N is odd. then the node 35 sits in the 
veiycemeiof the neighboihood. Oihervose. it sits slightly to one side. If two^ensional 
quincunx subsampling ^ perfonned, for even nmnbered levels, the local area is as shown 
in Fig. lOA. For odd numbered levels, the neighborhood is an NxN diamond, as shown in 
Fig. lOB. 

Referring back to Figs. 8 and 9, in step 82, the compaction tmit 28 selects a 
fust node 35 of the non-null pruned forest as the cunent "nucleus". In step 86 a zero- 
cicssing location for the cunent local area sunounding the cunent nucleus is computed If 
the location of the zero-crossing happens to be at the cunent nucleus (step 88), then the 
cunent nucleus is "zero-crossed" and there are no more iterations. Otherwise, the cunent 
nucleus is set to be at the location of the zero-crossing and the process is repeated for the new 
local area around the cunent nucleus. The iterations continue until the cunent nucleus is at 
the location of a zero-crossing or a maximum number of iterations (typically 5 to 10) has 
been reached. The last value for the zero-crossing, no matter how far from the original node, 
is saved (step 94) and the process repealed for each node of the non-null pnmed forest 

The equations for detennining the location of the zero-crossing are provided 
hereinbelow in the section labeled "Zero Crossing Equations". 

The remaining steps of Fig. 8 create "zero-crossed" uees which contain only 
those nodes which are zero-crossed or wUch have a zero-crossing within their local area. 
SpecificaUy, for each non-null pnmed tree, the compaction unit 28 scans the first resolution 
signal HI. Whenever (step 102) it fmds a "zero-crossed node" (i.e. one on which the zero- 
crossing occurs directly) from the list produced in step 90, it checks if the node is already a 
node in some zero-crossed tree (step 104) in which case, it stops the present detenninaiions 
and adds a temporary list STL to the zero-crossed tree to which the node belongs (step 1 10). 
If it is not already part of the zero-crossed uee, the unit 28 saves the node b the temporary 
list STL (step 106). 

Unit 28 then sets the cunent node to its parent node (step 108) and then has 
to determine if the parent node should be added to the zero-crossed tree. It repeats steps 104 
and 106 and, if the parent is not a root, the foDowing criteria are checked: (step 112) is the 
parent node zero-crossed, (step 114) does it have a zero-crossing in its local area based on 
the infonnatioh from the method of Rg. 9; or (step 116) are its parent and its child both zero- 
CTOSscd. In other words, step 116 indicates that the parent is a small gap in a branch of zero- 
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aossed nodes aad is a lelaxation criterion vis-a-vis the other criteria. 

If the paient node fulfills one of the criteria, it is added to the temporaiy list 
STL and the next pajent.is taken. The detennination is repeated imtfl a datapoint from the 
loot signal is leached (step 118) oi v^tn a node which is already part of a zero-aossed tiee 
is found. Otherwise, the detennination foi this node of the fust high frequency signal level 
HI is finished and the next node is taken. When the nodes of HI are exhausted, the process 
is repeated for the other high frequency signals, starting from level H2. 

The nodes are scanned in quanta of local area size. Thus, the first resolution 
level HI is divided into units the size of one local area. For a two-dimensional input image, 
the units are squares or diamonds for the odd levels from quincunx subsampling. During a 
first pass, the daiapoints in the same position in each local area, for example in the top left 
corner, are checked. In the next pass, the datopoints from the next position, for example, the 
next position in the top row, are checked. 

The datapobts are first checked that they are nodes of the non-null uee. If not, 

the next datapoint (i.e. from the next local area) is taken. 

The result of the operations of Fig. 8 are zero-crossed trees having therein only 
those daiapoints which are at or near transients. However, a human typically views a 
uansient as having a certain thickness. Thus, the zero-crossed trees do not contain enough 
information. They must be filled with those nodes of the non-nuU pruned trees which arc 
near the nodes of the zero-aossed trees. 

Therefore, in step 74 of Fig.'7, the compaction unit 28 adds data to the nodes 
of the zero-crossed trees. The operations performed in step 74 are listed, in How-chart 
format, in Rg. 11. 

The nodes of the new trees, caUed herein "ftilly pruned trees", are first set to 
be the nodes of the zero-crossed Uees (step 130). 

The compaction unit 28 then scans the nodes N of each non-null pruned tree 
to determine which nodes are either already part of the sinicture (step 132) or are part of the 
local area of at least one node of the conesponding zero-crossed tree (step 134). A list M 
is produced of the nodes of the zero-crossed trees whose local area the cuncnt non-null node 
N is in. 

If a node is neither part of the stnicture or part of the local area of a node, a 
new node is tested (step 136). If M has one member, m', the cunent node N is added to the 
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fully pnmcd uce at node m' (step 138). If M has many members O-c. the node N is part of 
many local aicas), the node m which is closest to the cunent node N is selected as m' (step 
140) and cunent node N is added to the fully pruned tree at node m* (step 138). Tht process 
is then repeated for all nodes in each non-null pruned tree. 

The process of Rg. 11 ensures that, if a node of a non-nvJI tree should be 
iBcluded in the fully pruned tree, it is included at only one node. Thus, the fully pruned tree 
contains the areas of high information content in a non-redundant manner. 

The convener 40 (Fig. 3) of the decompression unit 22 reproduces the wavelet 
transfonn data W* by placing the value of zero at every pixel not included in the list The 
inverse wavelet transformer 42 then reproduces the input sipal bom the uncompacted 
wavelet transform data W*. 



Nullifier Eoua lions; 



».(/^^>l5^*0J>J (5) 



where Ixl is the floor of X which is the maximum integer smaDer than X. 

Wjmj)^himW^{Uj\W^^ (6> 



Tl^y^-)..^ (8) 
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serve as 



iWil^j) otherwise 

Where RiUj) is a function of "imponance" which may be fixed for each level 
( /?(W>J«(0 ) 01 dependent (for example, the human visual response ftmcUon may 

ihe function R for image data). Q is a predetermined quantization step and W^^ 
is a predetermined maximum number of bins in the histogram. 
Zero Crossi tip; Fx|uations; 

The foUowing is noted about the zero crossing equations: a) they are written for a two- 
dimensional input signal, b) as the zero crossing is done within the level, we drop the level 
index, c) W(ij) denotes a detailed signal datapoint and d) the OperaUons are done for a w 
within a local area nucleus. 

lir>0 

ZZ-Ej'^^-'j^ ^^^^ 



(13) 
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(15) 



(16) 



(17) 



li will be apprecialed by peisons skilled in the art that the present invention 
is not limited to what has been particularly shown and described hereinabove. Rather the 
scope of the present invention is defined by the claims which foUow. 
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CLAIMS 



1. A system for compaction and leconsiniction of wavelet transform data of an 
input sipal, the wavelet Uansfonn data foiroed of a plurality of resolution levels, the system 
comprising: 

a compaction unit, operative on foiests of trees, each tree having branches 
which relate one daiapoint in a low resolution level with the entirety of its descendant 
daiapoints in higher resolution levels, for pruning branches from said trees and trees from said 
forest to remove aU but areas of high infoimaUon content in said input signal and for 

producing thereby pruned trees; and 

a convener for convening said pruned uees to reconstructed wavelet transform 

data. 

2, A system according to claim 1 wherein said compaction unit comprises a non- 
null pnming unit for removing branches of said uccs having daiapoints thereon whose wavelet 
transform values are below predetennined threshold levels and for producing thereby non-nuD 
pruned trees. 

3^ A system according to claim 2 and wherein said compaction unit additionally 

comprises a zero-crossing pruning unit for deiermining locations of transients in said input 
signal and for creating zero-crossed pnmed trees having therein only datapoints of said non- 
null pruned Uees which are at or near said locations of transients. 

4 A system according to claim 3 wherein said datapoints of said zero-crossing 

pnmed uees fonn nodes and wherein said compaction unit further comprises means for 
associating datapoints of said non-nuD pnined Uees with nodes of said zero-crossed pnmed 
uees, wherein each datapoint can be within a local area of many nodes but is associated with 
only one of them. 

5. A system according to claim 2 and including a nullifier, operative prior to said 

compaction unit, for nulling those datapoints of said wavelet uansform date which are below 
said predetennined threshold levels. 
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6. A sysiem according to claim 3 and wherein said zero-crossed pruning unit 

includes means for determining locations of zero-aossings in said wavelet transform data. 

5 7. A method of compaction and reconstruction of wavelet transform data of an 

input signal, the wavelet transform data formed of a plurality of resolution levels, the method 

comprising the steps of: 

pruning branches &om Uees and trees from forests to remove all but areas of 
high information content in said input signal and producing thereby pruned trees, wherein 
10 each uee has branches and each branch relates one datapoint in a low resolution level with 
the entirety of its descendant datapoints in higher resolution levels; and 

converting said pruned uees to reconstructed wavelet transform data. 

8. A method according to claim 7 wherein said step of pruning comprises the step 

15 of removing those branches of said uees having datapoints thereon whose wavelet uansform 
values arc below predetermined threshold levels and producing thereby non-null pnmed uees 

9 A method according to claim 8 and also comprising the step of determining 

locations of Uansients in said input signal and creating zero-crossed pruned trees having 
20 therein only datapoints of said non-null pruned uees which are at or near said locations of 
Uansients. 

10. A method according to claim 8 wherein said datapoints of said zero-crossed 
pruned uees form nodes and wherein said step of pruning further comprises the step of 

25 associating datapoints of said non-null pnmed uees with nodes of said zero-crossed pnmed 
uees. wherein each datapoint can be within a local area of many nodes but is associated with 
only one of them. 

1 1. A method according to claim 8 and including the step, operative prior to said 
30 compaction unit, of nulling those datapobts of said wavelet Uansform data which are below 

said predetermined threshold levels. 



20 



12. A method according to claim 8 and wherein said step of detennining includes 

the step of detennining locations of zero-aossings in said wavelet transfonn data. 
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